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|Let/r/^,  ^  be  k  independent  normal  populati 

f9Y, respectively,  and  a  common  known  variance  The  population 

*  -■  ?VA  t 

associated  with  the  largest/^  Is  called  the  best  population.  In  the  subset 
selection  approach, -we- Want  to  select  a  nonempty  subset  of  the  k  populations  so 
that  it  includes  the  best  population  with  a  minimum  guaranteed  probability. 

The  basic  idea  of  the  subset  selection  approach  is  that  the 
number  of  populations  to  be  selected  should  depend  upon  the  evidence  supplied 
by  the  data.  The  size  of  the  selected  subset  depends  on  the  sample  size  and 
the  confidence  level  P*  associated  with  the  claim  that  a  correct  selection 
(i.e.  selection  of  any  subset  that  contains  the  best)  occurs.  It  can  be  said 
that  subject  to  ^he  P*-value  and  the  sample  size,  we  cannot  make  finer  distinction 
among  the  populations  that  are  selected  in  seeking  the  best.  In  this  case,  one 
may  decide  to  use  the  selected  populations  In  equal  proportions  in  the  future. 

In  this  sense,  the  average  worth  of  the  selected  subset  is  given  by 


r. 


i  5*-l>  /< 


fLU  ?*■}>  k 

I.  INTRODUCTION 


ons  with  unknown  means 
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0.1) 


where  S  denotes  the  set  of  Indices  of  the  selected  populations,  and  ^  =  1  or 
0  according  as  tt..  is  or  is  not  included  in  the  selected  subset.  Our  interest 
is  to  estimate  H,  which  we  call  the  mean  of  the  selected  subset.  It  is  important 
to  note  that  M  is  a  random  variable. 

In  this  paper,  we  consider  the  subset  selection  rule  of  Gupta  [4],  [5], 
which  has  known  optimality  properties;  see  Gupta  and  Panchapakesan  [6],  and 
Berger  and  Gupta  [1].  Let  Y-j ,  Y2,...,  Yk  be  the  sample  means  based  on  n 
Independent  observations  from  each  population.  The  rule  R  of  Gupta  [5]  is: 


"Select  tt.  if  and  only  If  Y.  >_  max  Y.  -  dx/v^n  " 

l£jlk  J 

where  d  >  0  is  to  be  determined  such  that  the  probability  of  a  correct  selection 
(PCS)  is  P*.  This  value  of  d  is  shown  to  be  given  by 


(1.2)  f  *k-1  (t+d)  <f  (t)  dt  =  P*. 


where  (here  and  in  the  sequel)  $  and  <p  denote  the  standard  normal  cdf  and 
density  function  respectively. 

Our  present  Investigations  relate  to  only  the  case  of  k  =  2.  The 
following  notations  hold  for  the  entire  paper: 


(1.3)  *S 


*  2x2/n 

;  c  *  d/Y7  ; 

«  *  (ej-02)/o  ; 

“  Ve2 

;  0*  *  ; 

a  VY2 

;  Y*  =  Y^Yg  ; 

*1  *  I{Y>Co} 


lZ  *  !{Y<-co}  ; 


where  I.  denotes  the  Indicator  function  of  the  set  A. 
A 


Now,  for  k  =  2,  we  get 


c  =  ♦ 1 (P*) , 

(1.4)  ei+0? 

m  » 0 1 1 -j  +  02i2  +  -4-1  o-i1-i2)  • 

When  c  =  0,  the  rule  R  selects  the  population  that  yields  the  largest  sample 
mean  and  M  is  the  mean  of  the  selected  population.  Of  course,  in  this  case, 
the  minimum  PCS  cannot  be  guaranteed  for  P*  >  1/k  unless  additional  modifications 
are  made  in  the  formulation  of  the  selection  problem.  This  is  the  aspect  not 
considered  by  those  who  discussed  the  estimation  of  M  in  this  case;  these  are 
Sarkadi  [8]  and  Dahiya  [3]  for  k  *  2  and  known  x;  Hsieh  [7]  for  k  *  2  and 
unknown  x;  and  Cohen  and  Sackrowitz  [2]  for  k  ^  2  and  known  x. 

For  any  estimator  M  of  M,  the  bias  B(M)  -  E(T-M)  and  the  mean  squared 
error  MSE(M)  *  E(M-M)2.  It  can  be  shown  (Theorem  2.1)  that  no  unbiased  estimator 
of  M,  having  a  finite  variance,  exists.  In  Section  II  we  define  the  'natural' 
estimator  T  and  three  classes  of  estimators  T^x,  T2x  and  T3x,  obtained  by  making 
adjustments  for  the  bias  of  T.  The  biases  and  the  MSEs  of  Tlx,  T2x  and  T3x  are 
discussed  in  Sections  II  and  III  respectively.  Numerical  comparisons  of  the 
performances  of  these  estimators  are  made  in  Section  IV. 

II.  THE  ESTIMATORS  AND  THEIR  BIASES 
2 

Since  Y  *  N(e,o  ),  it  is  easy  to  see  that 

(2.1)  E (M)  »  9]{l-«(c-6)}  +  02{l-t(c+fi)>  + 

6  +0 

{$(c+6)+  *(c-6)  -1} 

=  J  [e*+  <5o{*(c+6)  -  <t>(c-6)}]. 


Theorem  2.1 .  No  unbiased  estimator  of  M,  with  finite  variance,  exists. 
Proof.  See  Appendix  A1.  D 

Let  us  now  consider  the  'natural'  estimator 

(max<Y1  ,Y2)  if  |Y|  <  co, 

Y* 

j  otherwise. 

Since  maxfY-pYg)  =  ^{Y*  +  |Y|},  we  obtain 

e(t)  =  f*  +  \  /  |y|  l  <f  {^)  dy  +  \  /  |y|  }  <r  (*~)  dy  • 

-00  Co 

By  changing  the  variable  of  Integration  by  setting  t  =  (y-e)/o  and  using 
(A2.1)  in  Appendix  2,  E(T)  simplifies  to 

(2.3)  E(T)  ■  |-  +  {«(c+6)  -  *(c-5)>  +  J  {q  (c+6)  +  q  (c-6 ) > 

=  E(M)  +  j  {tf  (c+6)  +  q  (c-6)}  . 

Thus  B(T)  ■  j{q(c+«)  +  q  ( C— 6 ) } . 

2.1  Estimators  T^ 

Since  the  bias  of  T  is  positive,  we  define  Tlx  by 
(2.4)  Tu  -  T  -  ^(T(c4)  ♦  t(c-J)) 
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X  >  0. 


6 


We  note  that  the  sum  of  the  first  two  terms  In  (2.7)  tends  to  E(M)  as  x  -*■  «. 
Also,  it  Is  seen  clearly  from  (A2.7)  that  the  last  term  in  (2.7)  is  unbiasedly 

estimated  by  y-  [q-  {X(^-  ♦  c)}  ♦  q  {k(~  -  c) }].  By  subtracting  this  unbiased 
estimator  from  U,  we  define 

(2.8)  T21  *  j  |  C*U(c*J))  -  ♦U(c-i))] 

-y  T  (x(c-X)J J. 

The  bias  of  T2a  Is 

(2.9)  B(T2a)  =  |  +  *  "  *  (6+c)  '  *  (6'c)j  ’ 

which  tends  to  zero  as  a  tends  to  Infinity.  Finally,  we  note  that  T2a 
corresponds  to  tc  of  Dahiya  [3]  when  c  *  0;  Dahiya's  c  corresponds  to  our  x. 

2.3  Estimators  T3a 

Let  us  first  consider  T3,  the  MLE  of  E(M) ;  this  Is  same  as  U  in  (2.6) 
with  X  =  1.  Thus,  from  (2.7)  with  X  *  1, 

(2.10)  E(Tj)  •  I*  4  |  r«(i±£)  -  4(^l)j 

2*7  L  *7  I 


Hence  the  bias  of  T3  Is 


(2.11)  B(T3)  =  |  +  '  $(<S+c)  "  $(6_c)] 

*-3-  U(~)  +  f  (^)1- 

2ft  L  ft  ft  J 


Noting  that  the  last  term  in  (2.11)  is  unbiasedly  estimated  by 


I  +  <F  (J  -  c)J  »  ^  define 

(2.12)  T3x  -  T3  -  x[v  j*(_|+ 


*(—  +  JL)  +  *(J( - £)  .  *(I  +  c)  -  *(J-  -c) 

i/?  i^b  /2b  ft 


♦|j*  (;♦*)♦, (I-oj]. 


The  bias  of  T3x  is 


(2.13)  B(T3  )  =  B(T3)  -  fir  (^)  +  if  (^)1 

JA  *  Zft  l  ft  ft  J 


-  4e  fv  j  *(-£  +  ~ )  +  *(■ 
2  L  (  SZc 


--  *  ~)  +  - -  *(j  *  C)  -  ♦<*  -  C) 

ft  fta  fta  ft  °  ) 


Now,  using  (A2.8)  to  evaluate  the  expectation  in  (2.13)  and  carrying  out  routine 
manipulations,  we  obtain 


^  *  •  v«.  *  *  O  *  ' 


(2.14)  B(T31)  ■  |  [(Hx) 


{(Si£)  +  ,(4S)|  .  x'.t(*S)  *  «t±% 

^  SZ )  (  /5  SS 


-  {*( 6+c)  +  3>(6-c)}j 


It  should  be  noted  that  is  slightly  different  from  the  estimator  of 
Dahiya  [3]  when  c  -  0  because  he  does  not  estimate  part  of  B(T3)  unbiasedly 
as  we  do  in  (2.12). 

III.  MEAN  SQUARED  ERRORS 

We  give  here  the  expressions  for  the  MSEs  of  T  and  T^.  For  and  T3a 
the  derivations  become  more  tedious.  For  the  numerical  comparisons  of  the  MSEs 
we  use  numerical  integration  as  will  be  explained  later  in  this  section. 


3.1  Mean  Squared  Error  of  T 


We  note  that  M  and  T  can  be  written  as 


(3.1) 


M  *  \  {e*  +  e(Ij-I2)}  , 
T  *  \  (Y*  +  Y(IrI2)>  . 


Now,  M  -  T  *  j  {(Y*-e*)  +  (Y-e ) ( 1^ -Ig) > .  Since  Y  and  Y*  are  Independent  and 
E(Y*-e*)  =  0,  it  is  easy  to  see  by  direct  evaluation  that 


(3.2)  MSE(T)  = 


t  +  t  /  tS(t)  dt  +  J  /  t*  tf  (t)  dt 

C-6  -• 

2 

=  {2  +  (c-6)  <p  (c-6)  +  (c+6)  (c+6)  +  4>( 6-c)  -  4>(6+c)}, 


using  (A2.3). 


3.2  Mean  Squared  Error  of  Tjx 
Y  Y 

Letting  V  =  <p  (c  +  -)  +  ?  (c  -  -),  we  have 

(3.3)  MSE(TU)  -  E[ (T-M)  -  %  V]2 

«  MSE(T)  +  (Xo2/4)  E (V2)  -  XoE[(T-M)V]. 

By  repeated  applications  of  (A2.6),  E(V2)  can  be  evaluated  in  a  straightforward 
manner  to  yield 

(3.4)  E(v2)  =  [<r  j\/f  («+c)|+  *  |y|"(5-c)|  +  ]. 

Again,  noting  that  V  (which  is  a  function  of  Y)  and  Y*  are  independent,  it  is 
easy  to  see  that 

(3.5)  E[(T-M)V]  ■  \  E[(Y-e)I1V]  -  \  E[(Y-e)IzV]. 

The  right-hand  side  of  (3.5)  can  be  written  as  a  sum  of  four  integrals 
each  of  which  is  either  of  the  form  in  (A2.4)  or  in  (A2.5).  Thus  we  get 
4  E[(T-M)V]  s  B,  +  B?  +  B,  +  B-,  where  # 
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Now  we  can  obtain  MSE  (T^x)  by  using  (3.2),  (3.4),  and  (3.6)  in  (3.3). 

Remark  3.1.  For  i  =  1,  2,  3,  T\ A  -  M  is  of  the  form  ^(Y*-e*)  +  g(Y), 

2  2 

where  g  is  some  known  function.  Hence,  MSE(Tix)  =  a  /4  +  E[g  (Y)];  the 
expectation  of  the  product  term  is  zero.  After  suitable  change  of  variable, 

E[g  (Y)]  can  be  numerically  evaluated  using  the  Gauss-Hermite  quadrature 
formula.  Our  results  in  the  next  section  were  obtained  by  using  the  20-point 
formula. 

IV.  COMPARISON  OF  THE  ESTIMATORS 

In  this  section  we  make  some  comparisons  of  the  performances  of  Tjx,  T^x 

and  T3x  based  on  the  values  of  their  biases  and  the  MSEs  computed  in  units 

of  a  and  o2,  respectively,  for  6=0  (0.5)  4  (1)  7;  c  =  1.645,  1.9600,  2.576; 

2 

X  =  0,  exp  {-c  /4> ,  1,  These  biases  and  MSEs  are  given  in  Tables  1A 

2 

through  1C  and  Tables  2A  through  2C  respectively.  For  convenience,  exp  {-c  /4> 
is  denoted  by  x  in  the  tables. 

V 

Remark  4.1.  The  choices  for  values  of  X  are  based  on  the  following 
considerations.  For  x  =  0,  Tlx  becomes  T.  The  value  x  =  1  corresponds  to 


11 


using  MLEs  in  constructing  the  estimators.  Further,  B(T^X)  decreases  in  x 
and  if  we  let  x^  to  be  the  value  for  which  B(T^X)  =  0,  Xfi  <  /2.  Also,  \& 
tends  to  0  or  •/!  exp  {-c  /4}  according  as  6  tends  to  infinity  or  zero. 


Remark  4.2.  The  c-values  chosen  here  correspond  to  P*  =  0.90,  0.95,  0.99, 
the  usual  values  of  interest  in  selection  problems.  The  value  of  c  =  0 
(P*  =  0.50)  is  not  of  interest  in  our  selection  problems.  This  is  the  case 


considered  by  Dahiya  [3].  As  pointed  out  earlier,  T^x  and  T2x  in  this  case 
coincide  with  M  and  t  of  Dahiya.  For  T,  ,  our  choices  for  values  of  x  are 

K  C  I  A. 

included  in  Dahiya 's  tables.  For  T2x»  X  =  1  is  the  only  common  choice.  However, 
we  do  not  report  our  values  in  this  case,  as  it  is  not  of  main  interest  here. 


Finally,  Dahiya  defines  a  hybrid  estimator  Hc  for  his  problem;  however,  Hc  is 
really  our  T.  But  his  c-values  are  chosen  arbitrarily  and  they  do  not  correspond 


to  P*-values  of  common  interest. 


Now,  considering  the  biases,  we  see  that  T2x  performs  better  than  Tlx  and 

T2x  for  small  values  of  6  (6  <_  1  for  P*  =  .90  and  6  £  .5  for  P*  =  .95,  .99). 

As  6  increases,  T^x  becomes  increasingly  bad  without  adjustment  for  bias  (i.e.  x  = 

however,  with  x  increasing  the  bias  of  T2x  is  very  much  reduced.  For  large 

6  (6  £  5  for  P*  =  .90,  6  >  6  for  P*  =  .95,  .99),  T^x  performs  better  than  T3x 

and  is  generally  better  than  T2x  as  well.  For  small  6  (6  <_  1  for  P*  =  .90,  .95) 

and$£  2  for  P*  =  .99)  T^x  performs  better  than  T3x  for  X  £  X£.  For  moderate  6, 

T,,  is  better  than  T,.  for  x  <  x„. 

3x  IX  —  c 

From  the  point  of  view  of  MSE,  an  overall  picture  emerges  as  follows: 


For  6  <_  1 ,  T2x  is  the  best.  For  1.5  £  6  £  3,  T^x  is  the  best.  When 

3.5  £  6  £  4,  T^x  with  small  x  or  T3x  with  large  x  is  the  best.  Finally,  for 

6  1  5,  T3x  with  x  away  from  zero  is  the  best. 


Table  1C.  Biases  of  T-jx(top  entry),  T2X(middle  entry), 
and  Tgx (bottom  entry)  expressed  in  units  of  a 
P*  =  .99;  c  =  2.576;  Ac  =  .2692 


\ 

6 

0 

Xc 

1 

.0145 

.0000 

-.0392 

0 

0 

0 

.0537 

.0332 

-.0225 

- 

.0249 

.0084 

-.0364 

0.5 

-.0045 

.0162 

.0096 

.0709 

.0479 

-.0144 

- 

.0579 

.0360 

-.0236 

1.0 

-.0287 

.0537 

.0347 

.1163 

.0874 

.0090 

- 

.1119 

.0828 

.0040 

1.5 

-.1057 

.0782 

.0603 

.1682 

.1333 

.0386 

- 

.1690 

.1338 

.0384 

-.2823 

.0411 

.0590 

.1896 

.1527 

.0527 

- 

.1989 

.0578 

.1609 

— 

2.5 

-.5871 

.0109 

-.0889 

.1519 

.0323 

.1197 

- 

.1823 

.1460 

.0474 

_ 

3.0 

-.9963 

-.2909 

-.0696 

- 

.0653 

.0438 

-.0145 

- 

.1302 

.0995 

.0162 

3.5 

-1.4389 

-.4978 

-.1383 

- 

-.0243 

-.0323 

-.0539 

- 

.0724 

.0495 

-.0126 

_ 

4.0 

-1.8456 

-.6442 

-.1595 

- 

-.0746 

-.0706 

-.0596 

- 

.0106 

.0018 

-.0219 

- 

5.0 

-2.4808 

-.0565 

-.7028 

-.0429 

-.0890 

-.0058 

.0006 

-.0015 

-.0070 

• 

-2.9991 

-.0148 

-.5980 

-.0060 

-.0223 

.0178 

* 

.0000 

-.0003 

-.0010 

7.0 

-3.5000 

-.0020 

-.4602 

.0010 

-.0031 

.0091 

• 
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Table  2A.  MSEs  of  T^(top  entry),  T2x (middle  entry), 

2 

and  T^x (bottom  entry)  expressed  in  units  of  a 
P*  =  .90;  c  =  1.645;  Ac  =  .7190 


A 

0 

A 

1 

<S 

c 

.3825 

.3605 

.3580 

.3606 

0.0 

.2500 

.2932 

.2925 

.2935 

.3195 

.3704 

.3962 

.4403 

.3841 

.3596 

.3564 

.3582 

0.5 

.2614 

.2917 

.2942 

.2997 

.3196 

.3726 

.3999 

.4466 

.3836 

.3566 

.3529 

.3545 

1.0 

.3288 

.2992 

.3100 

.3273 

.3256 

.3795 

.4085 

.4594 

.3778 

.3575 

.3566 

.3625 

1.5 

.5318 

.3418 

.3622 

.3916 

.3525 

.3952 

.4214 

.4699 

.3783 

.3764 

.3822 

.3973 

2.0 

.9402 

.4332 

.4570 

.4882 

.4107 

.4265 

.4434 

.4791 

.3983 

.4188 

.4320 

.4568 

2.5 

1.5628 

.5508 

.5652 

.5818 

.4878 

.4697 

.4732 

.4891 

.4337 

.4691 

.4865 

.5157 

3.0 

2.3499 

.6464 

.63  98 

.6310 

.5516 

.5061 

.4974 

.4941 

.4676 

.5051 

.5218 

.5485 

3.5 

3.2456 

.6857 

.6552 

.6244 

.5780 

.5198 

.5039 

.4877 

.4885 

.5184 

.5311 

.5509 

4.0 

4.2284 

.6710 

.6243 

.5854 

.5697 

.5128 

.4952 

.4738 

.4996 

.5100 

.5142 

.5206 

6.4994 

.5842 

.5412 

.5187 

.5247 

.4909 

.4789 

.4627 

.5000 

.5019 

.5027 

.5038 

6.0 

9.2500 

.5264 

.5071 

.5017 

.5042 

.4916 

.4868 

.4800 

.5000 

.5002 

.5003 

.5004 

FJFiK3&Si 

12.5002 

.5058 

.5007 

.5001 

.5004 

.4973 

.4961 

.4943 
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Table  2B.  MSEs  of  T^(top  entry),  Tg^middle  entry), 

2 

and  T3x(bottom  entry)  expressed  in  units  of  o 
P*  =  .95,  c  =  1.96,  Xc  =  .5413 


\ 


I 


I 


I 


\ 

6 

0 

xc 

1 

ft 

.3147 

.3090 

.3104 

.3166 

.2500 

.2802 

.2760 

.2755 

.2975 

.3204 

.3449 

.3709 

.3277 

.3191 

.3187 

.3240 

.2545 

.2806 

.2785 

.2798 

.3042 

.3316 

.3609 

.3920 

.3545 

.3396 

.3355 

.3387 

1.0 

.2902 

.2834 

.2882 

.2947 

.3201 

.3574 

.3975 

.4405 

.3726 

.3538 

.3480 

.3509 

1.5 

.4244 

.2992 

.3148 

.3286 

.3422 

.3846 

.4328 

.4859 

.3750 

.3606 

.3592 

.3663 

2.0 

.7500 

.3508 

.3744 

.3920 

.3779 

.4119 

.4565 

.5092 

.3779 

.3766 

.3851 

.4004 

2.5 

1.3282 

.4510 

.4691 

.4793 

.4358 

.4470 

.4743 

.5130 

.3983 

.4123 

.4315 

.4547 

3.0 

2.1404 

.5776 

.5697 

.5605 

.5052 

.4883 

.4915 

.5081 

.4337 

.4572 

.4820 

.5081 

3.5 

3.1082 

.6822 

.6337 

.6021 

.5585 

.5202 

.5026 

.4983 

.4676 

.4918 

.5149 

.5378 

4.1626 

.7296 

.6415 

.5966 

.5764 

.5298 

.5010 

.4836 

.4972 

.5087 

.5189 

.5285 

5.0 

6.4941 

.6797 

.5660 

.5061 

.5404 

.5061 

.4807 

.4608 

.5000 

.5026 

.5049 

.5070 

9.2500 

.5904 

.5140 

.5042 

.5087 

.4941 

.4824 

.4724 

.5000 

.5003 

.5006 

.5009 

■iOBB' 

12.5002 

.5359 

.5017 

.5002 

.5010 

.4969 

.4935 

.4905 
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Taking  the  bias  and  MSE  together  into  account,  the  relative  performances 
of  these  estimators  can  be  largely  summarized  as  follows:  For  6  ^  2,  T2a  is 
the  best.  For  2  <  6  <  5,  one  of  and  is  always  better  than  or  as  good 
as  T2x.  When  s  >_ 5,  T2a  is  never  the  best.  However,  for  6  >  2,  there  is  no 
clear  choice  between  TlA  and  T3a  because  the  bias  and  the  MSE  pull  each  in 
opposite  directions. 
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APPENDIX  1 

Proof  of  Theorem  2.1.  Since  Y*  and  Y  are  jointly  sufficient 

for  E(Mj,  we  restrict  our  attention  to  unbiased  estimators  that 

are  functions  of  Y*  and  Y.  Let  n(Y*,Y)  be  an  unbiased  estimator 

of  M  such  that  V(n]  <  ».  Let  Z  be  such  that  it  is  independent  of 

2  2 

Y*  and  Y,  and  Z  'v  N(6/b  ,  1/b  ),  where  b  >  0  is  known.  Applying 
the  Neyman  Factorization  Theorem  to  the  joint  distribution  of  Y*, 

Y  and  Z,  we  see  that  Y*  and  W  =  Z  +  Y/o  are  sufficient  for  0*  and 
6.  Thus  nj  =  E[n|(Y*,W)]  is  a  function  of  Y*  and  W  only.  Also, 
ECn^j  *  E(p)  and  V(p^)  <  V(.p)  by  Rao- Blackwell  Theorem.  Now,  let 

2 

n,  -  [—  +  -  Wf  {  4>  Cbc*  +  DW)  -  *  (be*  -  DW)}] 

L  L  1 

-[D  <p  (be*  +  DW}  +  D  q>  {be*  -  DW}] 

■  Ri  -  r2>  say, 

where  be*  =  c  /\  +  b^  and  D  ■  b^/  /l  +  b^ .  Using  (A2.7)  and 

(A2.8),  it  can  be  seen  that  E(n2)  *  E(M).  Since  ECn-^  -  n2)  =  °» 

it  follows  that  =  n 2  with  probability  1  by  the  completeness  of 

(Y*,w).  Thus  V(n2)  ■  v(n1)  <  V(n)  <  «. 

To  complete  the  proof,  we  obtain  a  contradiction  by  showing 

that  V(n2)  can  be  increased  indefinitely  by  letting  b  -*■  ®.  To  see 

this,  we  first  note  that  V(R^)  <  »  which  implies  that  V(R2)  <  «>. 

2 

It  is  also  easy  to  see  that  0  <  E(R2)  <  °°.  Hence,  E(R2)  <  00  which 

2  2 

implies  that  E(D  <p  {be*  +  DW}]  <  ®.  On  the  other  hand,  using  (A2.7) 

2 

and  the  fact  that  <p  (a)  *  <p(/2a )//!?,  we  can  see  that 


APPENDIX  2 


We  state  a  few  results  relating  to  normal  distribution.  We  omit  the 
proofs,  which  are  straightforward. 


(A2.1)  /  t  <r  (t)  dt  =  *  (a)  -  *  (b). 

a 


(A2.2)  / 


<f{oi(t+e)}  (t)  dt  =  — j—  (r  (-2$ — ^ 


(A2.3)  /  t  <f  (t)  dt  =  a  (a)  -  b  q-  (b)  +  $(b)  -  $( a), 

a 


(A2.4) 


00  r- 

/  t  <p  (t)  «r  (t+a)  dt  =  i  <r  (-£)  J ^  ( *^2a  +  —  )  +  —  *(/&  +  -S. 

a  ‘  ft  L  ft  ft  ft 


-b  ,  r 

(A2.5)  /  t  <r  (t)  *(t+a)  dt  -  -  I  ?(-2L)  L  (ft  b-  -*)+-*  *  (— 

*  ft  L  Si  ft  ft 


(A2.6)  /  <r  (t+a)  *  (t+B)  «f  (t)  dt  =  -1-  xf  . 

ftv  'V  J 


Let  Y  N(e,a^).  Then 


(A2.7) 

E 

c>|] 

.  >  „  jaliicU  , 

A+a2  (  A+a2  i 

(A2.8) 

E 

Ytjai-  +  Oil 

=  e  «  iaisislj  +  ™  v 

ja(6+c)  ) 

L  (  M 
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